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Abstract 

A dielectric columnar thin film (CTF), characterized macroscopically by a relative permittivity dyadic, 
was investigated theoretically with the assumption that, on the nanoscale, it is an assembly of parallel, 
identical, elongated ellipsoidal inclusions made of an isotropic dielectric material that has a different 
refractive index from the bulk material that was evaporated to fabricate the CTF. The inverse Bruggeman 
homogenization formalism was developed in order to estimate the refractive index of the deposited 
material, one of the two shape factors of the ellipsoidal inclusions, and the volume fraction occupied 
by the deposited material, from a knowledge of relative permittivity dyadic of the CTF. A modified 
Newton-Raphson technique was implemented to solve the inverse Bruggeman equations. Numerical 
studies revealed how the three nanoscale parameters of CTFs vary as functions of the vapour incidence 
angle. 

Keywords: Bruggeman homogenization formalism, Newton-Raphson technique, tantalum oxide, titanium 
oxide, zirconium oxide 

1 Introduction 

Columnar thin films (CTFs) are familiar structures within the optics literature, having been fabricated by 
physical vapour deposition methods for well over a century [1] . Their morphology is reminiscent of certain 
crystals, while their macroscopic optical properties are analogous to those of certain orthorhombic crystals. 
Furthermore, they are the precursors of the more complex sculptured thin films (STFs) [2]. 

The prospect of controlling the porosity and the columnar morphology of these thin films at the fabri- 
cation stage, in order to engineer their macroscopic optical responses, renders them attractive platforms for 
optically sensing chemical and biological species [3, 4, 5, 6, 7, 8, 9, 10, 11]. However, for intelligent design 
and deployment of such sensors, it is important to fully characterize the relationship between macroscopic 
constitutive properties on the one hand and the nanoscale morphology and composition on the other. 

There are significant impediments towards arriving at definitive relationships. One is the variability that 
exists due to differences in deposition conditions [12, 13]. For instance, the bulk material that is evaporated 
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may be quite different from the material that is actually deposited as a thin film. Therefore, while the 
dielectric properties of the bulk material is easily known prior to evaporation, the dielectric properties of 
the deposited material may well be different, depending on, whether the deposition occurred in an oxidizing 
or reducing atmosphere, whether trace amounts of water vapor were present, and the temperature. As an 
example, when Ti203 is evaporated, the deposited material has been shown by one research group to be 
either TiOi.g or TiOi.5, depending on whether the temperature is 25 °C or 250 °C [14]. Evaporation of 
different suboxides of titanium leads to the deposition of different TIOq films, in general, where the real 
number a varies with the nominal deposition conditions and even the deposition apparatus. Likewise, when 
Si02 is evaporated, the deposited material is some ill-defined but consistent mixture of Si and Si02 and is 
thus often classified as SIOq, a £ (1,2) [15, p. 164]. Furthermore, the delineation of nanoscale morphology 
is not an unambiguous task, as even a cursory glance at scanning-electron-microscope images of CTFs will 
confirm [1]. Direct determination of porosity or void volume- fraction through a gas-adsorption technique 
[16, 17, 18], although accurate, is very time-consuming. Therefore, porosity is usually measured indirectly 
through measurement of mass density, which has its own sources of inaccuracy [12]. 

Various researchers [1, 19, 20, 21, 22] have put forth nanoscale-to-macroscopic models for the relative 
permittivity dyadics of CTFs. Generally speaking, in these models the CTF is viewed as an assembly of 
parallel, identical, nanoscale inclusions of a certain shape dispersed in a certain homogeneous material. At 
optical and lower frequencies, these inclusions are electrically small and can therefore be homogenized into 
a macroscopically homogeneous material [23]. Apart from the shape of the inclusions, one must choose the 
porosity and the bulk dielectric properties of the deposited material and the material in the void regions 
(usually taken to be air) of the CTF. Such models require careful calibration against experiments [24]. 

Inversion of the forward homogenization procedure can provide nanoscale information about a CTF, 
which can be useful, for example, to predict what would happen if the CTF were to be infiltrated by some 
other material [25]. This thought motivated the work reported in this paper. 

Provided the components of the relative permittivity dyadic of a CTF are measured by suitable optical 
experiments [26, 27], an inverse homogenization procedure could yield the refractive index of the deposited 
material, the porosity of the CTF, and the shape of the inclusions, if certain reasonable assumptions are 
made. A demonstration based on the Bruggeman formalism [28] is presented in the following sections. In the 
notation adopted here, vectors are underlined whereas dyadics are double underlined. The unit Cartesian 
vectors are written as u^, Uy, and u^; the unit dyadic L = Ux2Ax + U.yUy + "MlzILz'^ tlie permittivity of free 
space is denoted by eo! the angular frequency is denoted by w; and i = \/— T. 



2 Homogenization model 

Let us consider a CTF grown on a planar substrate through the deposition of an evaporated bulk material. 
The planar substrate is taken to lie parallel to the z = plane, and the deposited material is assumed to be 
an isotropic dielectric material with refractive index Ug- At length scales far greater than the nanoscale, the 
CTF is effectively a continuum which may be characterized by the frequency-domain constitutive relation 
[1,2] 

D = eog^^p'E, (1) 

where 

^CTF ^ ^y(^) * {^<^^z^z + ^bU^Ux + ^cUylk^ ' i~^x) (2) 

is the relative permittivity dyadic of the CTF. The middle dyadic on the right side of Eq. (2) indicates the 
macroscopic orthorhombic symmetry of the CTF [1] . The orientation of the columns with respect to any xy 
plane is indicated via the inclination dyadic 

§,y (x) UyUy + {Ux2lx + u^Uz) cos X + (UzHx " ^xHz) sIh X , (3) 
where the column inclination angle is x G (0, 7r/2]. 
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Each column of the CTF may be regarded as a set of elongated ellipsoidal inclusions strung together 
end-to-end. All inclusions have the same orientation and shape. The latter is specified through the shape 
dyadic 

^U.n'lkn+ 7r Ikr + 76 U.b Ub^ (4) 

wherein the normal, tangential, and binormal basis vectors are specified in terms of the column inclination 
angle per 

Ikn = -ikx sinx + cosx ^ 

ILt ~ "iLx cos X + sin X > ■ (5) 

Ub = -y^y J 

Since the columnar morphology is highly aciculate, we have that the shape parameters 7& > 1 and 7t ^ 1. 
As increasing 7,- beyond 10 does not have significant effects for slender inclusions, we fixed 7r = 15 for 
definiteness. 

As the CTF is porous, we introduce / £ (0, 1) as the volume fraction occupied by the ellipsoidal inclusions 
representing the columns of the CTF. The void region is filled with air (or vacuum). Thus, the porosity of 
the CTF equals 1 - /. 



3 Forward and inverse homogenization 

The nanoscale parameters {us, /, 7;,} may be related to the eigenvalues {ta, Cc} of ^Qj^p ^'^'•^ one of several 
homogenization formalisms, including the Maxwell Garnett formalism [22], the Bragg-Pippard formalism 
[1] , and the Bruggeman formalism [2] . We implement here the last-named formalism which has been widely 
used in optics [29], because it treats the region occupied by the deposited material and the void region 
symmetrically, unlike the other two formalisms. 
Let us introduce the dyadic 

^ = /a, + (!-/)%' (6) 

which is the volume- fraction- weighted sum of the two polarizability density dyadics [2, 30] 
g, = eo [nl L - (ea u,^u^, + u,^.u^ + UyUy)] 

• |l + iweo^^ • [nlL~ {eaUz2l^ + <^bUxUx + <^cUyUy)]^ (7) 

and 



Co [L - {eaUz^z +^bUxikx + ^cHyUy)] 

• [L+i^<^0Rj: ' [L~ i<^aUzUz+<^bUxUx + <^cUyUy)]^ ■ (8) 

Herein, the depolarization dyadics 



1 9 /••^/2 

D = — - / d(p de 



ILxlLx + Sin^ I COS^ "ikzikz + '^'"a'^ ILyik^ 



sine (9) 

^+sin2 0(e,cos2 0-f-ee^^ 



and 



1 



= — — -I d<j) I dO 



it/2 j--iT/2 



uvea TT J^^Q Je^Q 

cos^ 9 u^u„ + sin^ 9 (cos^ (h u^u, + sit? 6 u„u,,) 

sin 9 =^=^ j^^T — -^=^ , (10) 

Cb cos^ 9 + sin 9 [Ea cos^ (/) + ec sin -'-^ 
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are straightforwardly evaluated by numerical means. According to the Bruggeman homogenization formal- 
ism, the three parameters {us, f,"fb} satisfy the three nonlinear equations 

bi{ns,,f,'yb) =0, {£^x,y,z), (11) 

where bi are the three nonzero components of the diagonal dyadic b ; i.e., 

k^bxU^u^ + byUyUy + bzU^u^. (12) 

Usually, the process of homogenization is applied in a forward sense, to provide a nanoscopic-to-continuum 
model. Thereby, the relative permittivity parameters {eq, eb, Cc} may be estimated from a knowledge of the 
nanoscale parameters {ris, /, 7t,}. However, the nanoscale parameters of CTFs are generally unknown whereas 
{ea,ef,,ec} may be measured. In order to estimate {n^,/, 7b} from a knowledge of {ea,eb,ec}, the inverse 
homogenization process is needed. Formal expressions of the inverse Bruggeman formalism are available [31], 
but in certain cases these formal expressions may be ill-defined [32]. In practice, it is more convenient to 
implement a direct numerical method to compute {'t-s,/, 7b}, as described in the next section. 



4 Numerical implementation 



Solutions to Eqs. (11) may be computed using a modified Newton-Raphson technique [33, 34]. In the 
recursive scheme implemented here, the estimated solutions at step fc -I- 1, namely |ni'^~''^\ f^''^^\ Tb'""''^^}' 

are derived from those at step fc, namely |ri,i'^-', Z^*^', 7^''''|, via 

bM'\f^'\ii'^) 



f(fe+l) 



df 



% 



(k+1) _ (k) 



% 



(fe)^ 



(13) 



^6.(ni'^-+i\/('=+i),7f 

In order for the scheme (13) to converge, it is crucial that the initial estimate |rii°\ /'°'', 7^''''| be 
sufficiently close to the true solution. A suitable initial estimate may be found by exploiting the forward 
Bruggeman formalism, as follows. 

Let ea,b,c denote estimates of the CTF permittivity parameters ta,b,c, computed using the forward Brugge- 
man formalism for physically reasonable ranges of the parameters n^, / and 7b, namely e (n^,n^), 
./ e (/^, n and 7b e {it.l'i!) ■ Then: 

(i) Fix Us = {n^ + n^) /2 and 7b = {-/t + lb) /2- For aU values of / e {f^J^), identify the value /* 
for which the quantity 

A = {ea - iaf + (eb - hf + (cc - icf (14) 

is minimized. 

(u) Fix f ^ f* and 7b = (7^ + 7^) /2. For aU values of e (nf", ), identify the value n* for which A 
is minimized. 

(iii) Fix / = /* and Ug 

The steps (i)-(iii) are repeated, using n* and 7^ as the fixed values of and 7b in step (i), and 7^ as the fixed 
value of 7b in step (ii), until A becomes sufficiently small. In our numerical experiments, we found that when 
A < 0.01, the values of n*, /* and 7^ provide suitable initial estimates for the modified Newton-Raphson 
scheme (13). 



I*. For all values of 7b S (7^ , 7^) , identify the value 7^ for which A is minimized. 
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5 Numerical results 



We considered CTFs made from three different materials: the oxides of tantalmn, titanimir and zirconium. 
Experimental studies [26] have revealed that the permittivity parameters for these CTFs may be expressed 



as 



(15) 



Ca = {naO + nal V + na2 V^) 
Cfc = ("'fcO + mi V + 7162 v'^Y 

Cc = (7^cO + nci f + nc2 v'^Y 
V = 2x«/7r 

wherein the vapor incidence angle Xv G (0i''''/2] is related to the column inclination angle by the coefficient 



TO, per 



tan X = TO tan Xv ■ 



(16) 



Table 1 contains values of the ten coefficients riao to to of CTFs of the three different materials. Although 
the bulk refractive indexes of all three oxides are quite close to each other, the coefficients Uao to to of the 
three types of CTFs are quite different, as indeed are also their constitutive parameters ea.b,c [35]. These 
differences arise, in significant measure, due to the dependence of the growth dynamics of a CTF on the 
evaporated bulk material [36, 37]. 



Table 1: Coefficients appearing in Eqs. (15), obtained from the experimental findings of Hodgkinson et al. 
[26] on CTFs, when the free-space wavelength is 633 nm. 



material 




nal 






nbi 


"62 


tantalum 
oxide 


1.1961 


1.5439 


-0.7719 


1.4600 


1.0400 


-0.5200 


titanium 
oxide 


1.0443 


2.7394 


-1.3697 


1.6765 


1.5649 


-0.7825 


zirconium 
oxide 


1.2394 


1.2912 


-0.6456 


1.4676 


0.9428 


-0.4714 



material 








to 


tantalum 
oxide 


1.3532 


1.2296 


-0.6148 


3.1056 


titanium 
oxide 


1.3586 


2.1109 


-1.0554 


2.8818 


zirconium 
oxide 


1.3861 


0.9979 


-0.4990 


3.5587 



The nanoscale parameters {n^,/, 7b} were estimated for the three CTFs using the modified Newton- 

Raphson technique (13), with initial guesses |rii°'', Z^^-*, 7b°''| deduced by scanning the solution space of 

{~ea, h, ed with = 1, = 4, = 0.2, = 0.9, = 0.5 and 7,^' = 3. 

The computed nanoscale parameters n^, /, and 7b, respectively, are plotted in Figs. 1-3, against Xv G 
(12°, 90°) for the CTFs fabricated by evaporating any one of the three bulk materials. The plots in Fig. 1 
show that Ug for CTFs made by evaporating any of the three bulk materials to be largely insensitive to Xv- 
In contrast, the volume fractions / displayed in Fig. 2 for all three materials increase rapidly as Xv increases, 
in general accord with the observation that mass density of a CTF varies as (1 + sinxw)"^ sinx^ [38]. The 
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Figure 1: The quantity Us plotted against Xv (in degree) for CTFs made from evaporating titanium oxide 
(red, solid curve), tantalum oxide (green, dashed curve) and zirconium oxide (blue, broken dashed curve), 
as computed using the inverse Bruggeman formalism. 
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Figure 2: As Fig. 1 except that the quantity plotted against Xv is /. 

shape parameters ji, displayed in Fig. 3 for CTFs of all three evaporated bulk materials decrease rapidly 
towards unity as Xv increases. This is in accord with the observation that CTFs deposited with Xv ~ 90° 
are macroscopically uniaxial rather than biaxial [1]. 

6 Concluding remarks 

In order to exploit the considerable potential that CTFs possess for widespread applications such as optical 
sensors of analytcs, it is vital that they be reliably characterized at the nanoscale. Our theoretical and 
numerical study has demonstrated that the inverse Bruggeman homogenization formalism provides a prac- 
ticable means for this characterization, in terms of three nanoscale parameters. Thus, a key step towards 
the intelligent design and development of CTF-based (and other STF-based [39]) optical sensors has been 
taken. 
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Figure 3: As Fig. 1 except that the quantity plotted against Xv is 7&- 
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